Thermal van der Waals Interaction between Graphene Layers 
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The van de Waals interaction between two graphene sheets is studied at finite temperatures. 
Graphene's thermal length (^t = hv/ksT) controls the force versus distance (z) as a crossover from 
the zero temperature results for jz <C ^t, to a linear-in-temperature, universal regime for 2 ^ ^t- 
The large separation regime is shown to be a consequence of the classical behavior of graphene's 
plasmons at finite temperature. Retardation effects are largely irrelevant, both in the zero and finite 
temperature regimes. Thermal effects should be noticeable in the van de Waals interaction already 
for distances of tens of nanometers at room temperature. 
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I. INTRODUCTION 

Graphene, the single layer honeycomb lattice of car- 
bon atoms that forms graphite, has been realized ex- 
perimentally in recent timesii^. Its electronic properties 
characterized by a linear dispersion around Fermi points, 
fixed by charge neutrality {massless Dirac fermions with 
velocity v ~ 10^ "m/s), have long attracted theoretical 
interestr^.. But it is the present experimental accessibil- 
ity, including Fermi level tuning by gate voltages, what 
has unleashed an explosion of activity, fueled in part by 
the prospects of tailoring its electronic (and perhaps mag- 
netic) properties in the nanoscale^. 

More traditional areas like the van der Waals (vdW) 
interaction have also benefited from the present inter- 
est. Although graphite is often characterized as a vdW 
stack of graphene layers, fundamental aspects such as 
the asymptotic behavior of the vdW interaction between 
two graphene layers, have been unveiled only recently 
by Dobson et al.^ , and shown not to conform to the 
naive sum of R^^ contribution^. Taking as reference 
the progress in accurate measurements of vdW interac- 
tions in general^, the expected increase in availability of 
graphen^ii, and its unique conceptual place as neither 
a metal nor a dielectrioS., the study of graphene's vdW 
interaction seems worth of further consideration^. 

In this paper I consider the vdW interactions between 
two graphene layers at finite temperature (T). Graphene, 
being a critical system at zero T, lacks any characteristic 
length scale'^. Temperature provides such scale, the ther- 
mal length: — hv /ksT . We will show that the thermal 
length controls the vdW interaction between planes in 
the form of a crossover. For separations between the two 
layers [z) smaller than the thermal length, z < ^t, the 
zero-T result (/) for the forced prevails, / cx 1/z^. But 
for separations larger than the thermal length, z > ^t, 
the force crosses over to a linear-in-T^'^ material pa- 
rameters independent, universal regime, / oc T/ z^, that 
constitutes the genuine asymptotic large-separation in- 
teraction between two graphene sheets at finite T. 

The linear-in-T regime will be shown to reflect the clas- 
sical nature of graphene's low lying excitations at finite 



T: plasmons. As shown by Vafek-ii, these plasmons are 
the charge fiuctuations of thermally generated carriers 
(electrons and holes). Therefore, they are present only 
at finite T and with energy scale tied to T, so that long- 
wavelength plasmons always behave classically. As such, 
this thermal regime will be shown to be present event 
for the instantaneous (non-retarded) Coulomb interac- 
tion. This should be contrasted with the usual linear- 
in-T, thermal limit of the vdW interactions between any 
material^^ that sets in for distances larger than the 
thermal length of the field, Xt = hc/ksT . The explicit 
appearance of the light velocity c in this generic case, 
is a manifestation of the classical population of the rel- 
evant electromagnetic modesi^. But not in graphene, 
where the existence of this regime even without retar- 
dation c — > cxD, and with the role of c taken by v in 
setting the range, shows it to be a consequence of the 
classical dynamics of matter. As a corollary, the inclu- 
sion of field's dynamics (retardation) will be proven to 
be largely irrelevant both at zero and finite T. This is 
another aspect where graphene separates from ordinary 
dielectrics and metals, where retardation always matters 
for large enough separations. The relevance of the ther- 
mal length (^t) (as opposed to At) places graphene in 
an unique position for the experimental observation of 
thermal effects in the vdW interaction. For instance, at 
room temperature, the thermal regime should begin to 
be observable for distances -z ^ ~ 26 nm. In contrast, 
the onset of thermal effects linked to the classicality of 
the field (the situation for good metal and dielectrics, see 
section |V| would require much larger distances at room 
temperature: z > At ~ 300 ^t- 

Concerning the experimental situation, a word of cau- 
tion is required. Throughout this paper, graphene is 
modeled as the usual set of Dirac fermions, known to pro- 
vide an excellent description of the low energy physics^ 
(even up to ~ 2 eV in light absorption experiments^^). 
Nevertheless, the vdW calculation is incomplete for real 
graphene, where more electrons and bands contribute to 
the vdW force, particularly at short distances. These 
additional terms are estimated from published ab initio 
resultsi^ in appendix [Bl and shown not to compete with 
the low energy contribution already for distances of the 



order of nanometers. 

The paper is organized as follows. In section [TTl the 
formalism is presented and the results^ at zero T recov- 
ered. Finite T is considered in section lllll the force versus 
scaled distance (z/^t) distance is computed numerically 
and its asymptotic behavior explained in terms of clas- 
sical plasmons. Full retardation is included in section 
IIVI where its contribution is shown to be quantitatively 
irrelevant for graphene's parameters. Section |Vl summa- 
rizes graphene's vdW results, contrasting them with the 
known behavior of ideal dielectrics and metals. Appendix 
|X] calculates graphene's transverse-response contribution 
to the vdW force, to show that it never competes with the 
the longitudinal response considered in the body of the 
paper. Appendix [B] addresses real graphene, estimating 
the contributions to the vdW force beyond that of Dirac 
fermions considered here. 



II. FORMALISM. ZERO TEMPERATURE 
RESULTS 
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FIG. 1: Template diagram for the computation of the vdW 
force, as in Eqs. 1171 and lAll Bubbles: isolated graphene's 
response (charge or current). Thin wavy lines: Coulomb inter- 
action (or photon propagator). Thick wavy lines: multiple- 
scattering-corrected interaction (or photon propagator) be- 
tween graphene's layers as in Eq. |3] Zig-zag line represents 
force (/c): distance derivative of the interaction line (or pho- 
ton propagator). 



Let US first present our non-retarded formalism recov- 
ering the zero-T resuHi^. Consider two graphene layers 
perpendicular to the z-axis and separated by a distance 
z. Ignoring (for the moment) retardation effects, the mu- 
tual force per area can be written as: 



/ = 



_fq_ 

(27r)2 



(1) 



with Coulomb coupling between density fluctuations pq 

and p*^^ given by Vc{q, z) = exp(-g|z|)/(2eog) (ele- 
mentary charge e and vacuum permittivity eo, SI units), 
and Coulomb force fc{q, z) = —dzVc{q, z). 

Evaluating the thermal average to all orders in the mu- 
tual interaction, we can write: 



/ = 



1 



(fq 



fcjq, z) 



i^n,z) xf'jiqjujn), (2) 



with = ksT, Matsubara frequencies fnon = 2TmkBT, 
and the multiple-scattering-corrected interaction be- 
tween planes given by: 



Wc{q,uj,z) 



Vc{q,z) 



1 - Vciq,z)^x'ppi<l,'^)x'ppi<l,^) 



, (3) 



where x'pp — Xp^p) = Xpp is the charge-charge Green's 
function of an isolated graphene layer, which can be writ- 
ten as: 



Xpp(9,w) 



Xpp [Q, w) 



1 - Vc{q,z^ 0)x%\q,^^) 



(4) 



with xlyj (q, '-^) as the polarization of an isolated graphene 
{proper polarization in diagrammatic sense^^). Before 



proceeding to the evaluation of /, let us remark that 
this formula is entirely equivalent to the (non-retarded 
version of) Lifshitz treatmenfc!'^*', as can be seen by eval- 
uating the field's stress tensor— in the presence of the 
(here non-local) material's response. If we knew the ex- 
act polarization of a single graphene (including its crucial 
q-dependence) , the only remaining approximation in ([2]) 
(and in Lifshitz's approach) would amount to the ne- 
glect of proper polarization diagrams connecting both 
planesi^: local field corrections to the dielectric response, 
safely ignored for large separations. The diagrammatic 
description of the vdW calculation is presented in Fig. [T] 
Throughout this paper, we will take as the proper po- 
larization the non-interacting value, what amounts to the 
standard RPA for Eq. (|4]). It is given by: 



X^J = N 



with N=2x2 fermion species, /^^ — ^ 
is the Fermi factor and E^ — ahvk. If we use the zero-T 
valueilii^ for x'"^ : 

X^"Hq,i^) 



N 



16hv y/q2 - uj2/v^ 



(6) 



in Eq. f!]), with the Matsubara sum becoming an integral 
iuJn — > irj, the resulting expression for / is: 



1 

2^ 



q dq 



drj ■ 

Jo 

exp(-2g|z|) 



16 
Na 



exp(-2g|z|) 



(7) 
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where a = e'^/{2eohv) is a dimensionless measure of the 
effect of interactions in graphene, with value a ~ 13.6. 
Expression [7] can be shown to be entirely equivalent to 
the treatment of Dobson et al^, leading to the follow- 
ing quantitative value for the force per area between 
graphene layers: 



A 



A - 0.40 eV A. 



(8) 



Notice that other choices for the proper polarization com- 
plying with the scaling Y^'^-' oc qf{Lu/vq), such as the ex- 
citonic response of ref. [l^, would produce the same z"^ 
power law, although with different prefactor. 



0- 



>" -0.2 - 



-0.4- 




III. FINITE TEMPERATURE RESULTS 

Let's consider now a finite temperature still for non- 
retarded interactions. Although no simple analytical 
expression is known for x^*^^ a,t finite T (see, though, 
ref.[20[). its scaling behavior— is best described measuring 
lengths in units of = hv/ksT, and energy in terms of 
ksT. Indeed, for the force calculation, matter and field 
appear in the dimensionless combination: 

v,{q,z)x'-"Hq,^)^a eM-<lM)xiq(T,hLu/kBT), (9) 

where xiq^T, fi^/ksT) is a dimensionless function. It is 
clear that the force, Eq. will depend on distance and 
temperature only through the combination zjS^r-, with 
the following scaling form: 



fiz,T) = ^fiz/^T) 

ST 



(10) 



Therefore, knowledge of the dimensionless function 
/{z/^^t) provides all information for the vdW interaction 
at finite T and arbitrary distances in the scaling regime. 
We have evaluated numerically the force (Eqs. (O, Q 
and ([U) with results plotted in Fig. [5J As expected, 
graphene's thermal length marks a crossover between 
two regimes: the zero-T limit ([5]) for z/^t ^ 1 previously 
analyzed?-., and the genuine large-distance regime at finite 
temperature for z/^t 3> 1, that we now consider. 

As discussed by VafekU., the most important feature 
of the charge-charge response ([3]) at finite T, is the emer- 
gence of plasmons. These appear as the zeros of the de- 
nominator of Eq. ^ and, in the long wavelength limit, 
the plasmon frequency is given b y^^i^^ : 



hx!p{q) 
keT 



In 2 
"2^ 



N a q ^T, <? f T ^ 1 



(11) 



Although plasmons possess an imaginary part 7(9), 
meaning that they decay into the electron-hole 
continuurftii, they become very long-lived excitations for 
long wavelengths: j{q)/ujp{q) — > for ^ 0. 

Being the density fluctuations of thermally excited car- 
riers, plasmons owe their existence and energy scale to 



FIG. 2: van der Waals (scaled) force per area / between 
graphene planes versus distance z in units of the thermal 
length ^T. Continuous line: numerical result. Dashed line: 
large distance {z 3> ^t) limit, Eq. (|14p . Dashed-dotted line: 
zero-temperature limit {z <^ ^t), (Eq. IS}. Inset: enlarged 
view of the crossover region {z ~ ^t). 



temperature. In this respect, they differ from plasmons of 
ordinary 2d metals, already present and contributing^;^ 
to the vdW force at zero T. Furthermore, the spectral 
power of the charge-charge response for hvq < huj < ksT 
and q^T ^ 1, is dominated by the plasmon mode. There- 
fore, a single-pole approach for the response suffices for 
the plasmon contribution to x ^-t finite T, with the ex- 
plicit form: 



X(9,^) = 



1 



ujpiq)'^ 



Vciq, z = 0) ^2 - LUp{qy 



(12) 



This response is valid for g ^ 1 and fku < ksT. Its 
use in Eq. ^ provides the wanted large-distance behav- 
ior of the vdW force at finite T. The evaluation is best 
performed trading Matsubara sums for real frequency in- 
tegration in Eq. 



1 

'2^ 



q dq 



dcu coth(/3Scj/2)- 
exp{-2q\z\) 



cxp(-2q|z|)' 



leading to the following central result: 



/ 



C(3) keT 



(13) 



(14) 



with the Riemann's zeta bmction, C(3) = 1.2020... . 
Eq. is the large-distance asymptotic behavior for the 
force per area of two graphene sheets at finite tempera- 
ture. The numerical solution does indeed merge with this 
analytical limit for z ^ ^t, as seen in Fig. [2l 

The result of Eq. is truly remarkable: all ma- 

terial and electrical parameters have disappeared, leav- 
ing the temperature as the only surviving energy scale. 



4 



As remarked in the introduction, an identical formula 
describes the force between two metallic plates at finite 
temperature, for distances larger than the thermal length 
of the electromagnetic field^°, z ^ At = ftc/ZcsT. This 
is the limit where the thermal population of the relevant 
electromagnetic modes becomes classical. But, in spite 
of the similarity, we cannot make an obvious connection 
with our result: our treatment has been obtained for the 
instantaneous, non-retarded Coulomb interaction, there- 
fore there is no field dynamics, no field modes, and the 
issue of classicality for the field is out of place. Setting 
c = oo in Ar renders meaningless the would-be range for 
that classical limit. Yet, our regime of Eq. p4)) for the 
instantaneous interaction appears for z > = hv/ksT. 

Nevertheless, the fact that v takes the role of c in set- 
ting the range for our non-retarded calculation prompts 
for the existence of a classical interpretation, but now 
for the only dynamical entity so far considered: matter. 
Plasmons, by the very fact that their existence and scale 
are tied to temperature, behave classically at long wave- 
lengths: 



0, 9eT«l, 



(15) 



and this suggests that there must be more transparent 
ways of getting such a simple result as Eq. ([M]) . As 
reassurance that our reasoning is well founded, we will 
now recover Eq. p4p invoking only elementary classical 
concepts. Let's consider graphene's charge fluctuations 
as classical objects at temperature T. The classical limit 
means that we can ignore kinetic energies and rely only 
on the potential (electrostatic) energy to account for the 
thermal population of these fluctuations. This electro- 
static energy is: 

= z)p^^)p^^l + U{q, 0)(p«p« + p^^^p%) 



(16) 



q (T— ± 



where we have diagonalized the quadratic form with 
the normal modes: p^^ — {l/V^){pq^ ± Pq''), with 
v±iq, z) = Vc{q, 0) ± Vc{q, z). The equipartition theorem 
allows us to write the thermal population of modes as 
< Pq^V'^'' >— ksT / v±{q, z). Expressing Pq'^^ in terms 

of Pq^'', the thermal average of Eq. ([T]) can be obtained 
with the result of Eq. HH). This fully supports our in- 
terpretation that it is the classical population of thermal 
plasmons what leads to the vdW force. 

Preparing for the discussion of retardation in the next 
section, it is worth noticing that the expression of Eq. [TJ] 
amounts to selecting the iojn = term in the frequency 
sum of Eq. [51 Indeed, this is expected for a classical limit, 
but notice that the reason for such behavior is entirely 
due to the matter response x and, again, has nothing 
to do with the interacting fleld which, considered so far 
instantaneous, therefore shows no frequency dependence. 



IV. RETARDATION 

Now we address the issue of the electromagnetic field 
dynamics, to show that the inclusion of retardation 
hardly affects the previous results. Following ref. [l^, re- 
tardation effects are best handled in a gauge where only 
the vector potential A exists. The coupling matter-field 
is of the form cx j- A, where the current lays in graphene's 
planes, and can be decomposed into (in-plane) longitu- 
dinal and transverse components that are not mixed by 
the photon field. Let's consider the longitudinal current 
responsible for charge fluctuations. It is straightforward 
to show that retardation can be included in the previous 
formalism with the following correspondences in Eqs. ^ 
and dSI): 



Vc{q,z) 
fc{q,z) 



A(g,w,2:) 
-dyVi{q,ijj,z) 



(17) 



where Vi{q,uj, z) is the (part of the) photon propagator 
that couples to in-plane longitudinal currents, with ex- 
pression: 



T^i{q,^,z) 



e^q' exp(— q'|z|) 
2eo cj2 ' 



(18) 



and q' ~ \J (f' ^ io^ I c? . Xjiji is the longitudinal current- 
current response, related by particle conservation to the 
charge-charge response by q^Xjiji (9i ^) = ^'^Xppilj It 
can be checked that setting the light velocity c — ^ oo 
in the above expressions, the non-retarded expression for 
the force is recovered. Again, as in the non-retarded case, 
this formalism for the force can be shown to be exactly 
equivalent to Lifshitz's when applied to the longitudinal 
response. 

Let's consider the results for zero temperature first. 
Carrying out the prescription of Eq. ([T7]) and trading q 
for g', the following expression is obtained for the vdW 
force due to longitudinal currents, including retardation: 



f = -^ q"dq'- dv 
27r Jo Jo 

exp{-2q'\z\) 



(1 + ^V^ + Vy{veffq'r ? " exp(-2g'|z|) 



(19) 



where v~^j = v^-^ — c^^. Comparing Eqs. 1191 and [71 one 
sees that only two formal changes appear with respect 
to the zero-T, non-retarded calculation. First, there is 
a renormalization of graphene's velocity in the square 
root of Eq. ([6]), u — s- v j \J\ — f^/c^, quantitatively irrel- 
evant. Second, the integration over imaginary frequency 
acquires an upper limit, whose physical interpretation 
corresponds to the removal of the electromagnetic field 
modes that, for each space scale, are slower than mat- 
ter. This is the dominant effect of retardation but, the 
integrand decaying as ~ ?7~^, it amounts to a meager 
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Dielectric Graphene 




q q 

FIG. 3: Left panel: light-cone (u = cq) versus dynamics of a 
dielectric of typical frequency ujq, illustrating the importance 
of retardation for distances beyond Zret ~ iTst- Right panel: 
light-cone versus graphene's typical dynamics (uj = vq, not to 
scale) illustrating the absence of a characteristic distance for 
retardation. 

v/c ~ (300)""'^ fractional reduction of the prefactor A in 
Eq. ([5]) without altering the power law. 

The irrelevance of retardation in graphene at zero 
T contrasts with the situation for a regular dielectric, 
where retardation always matters beyond some distance 
Zret ~ Iret^ with q^et ^ ^o/c, whcrc tujjo IS a typical en- 
ergy scale (say the gap). This is schematically depicted 
in the left panel of Fig. [S] In graphene (right panel 
of Fig. on the contrary, both matter and field are 
scale-invariant (critical) systems (with dynamical critical 
exponent 1), this implies that the ratio {c/v ^ 300) of 
their relative dynamics remains the same at every length 
scale (separation between planes). Therefore, the irrel- 
evance of retardation effects in graphene at zero T is 
both qualitative and quantitative. Qualitative because, 
at least within our RPA treatment, the power law for 
the vdW force of Eq. ([5]) would remain the same for 
arbitrary values of graphene's velocity w, although with 
a changed prefactor. For graphene, this irrelevance is 
also quantitative, because the prefactor barely changes: 
5 A/ A - 1/300. 

Now we show that retardation at finite T also lets un- 
affected Eq. (ITil) as the correct large distance behavior. 
The spatial dependence of the photon propagator ([18]) 
makes short-ranged the contribution from Matsubara fre- 
quencies other than n = 0. This effect begins to matter 
for distances 2 > Ay = hc/ksT, and is present in any ma- 
terial as it corresponds to the above mentioned classical 
limit of the thermal population of electromagnetic modes. 
But for graphene, restricting to n = adds nothing to 
the non-retarded result for Eq. p^ . Indeed, such result 
is equivalent to selecting n = in the matter response, al- 
though in that case this restriction was forced upon us by 
the classical behavior of matter's plasmons while the field 
remained instantaneous. In other words, for the vdW in- 



teraction in graphene, there is no difference between clas- 
sical matter + instantaneous field and classical field. We 
have computed the force with the numerically evaluated, 
finitc-T response of Eq. ([5]) and the retarded interaction, 
with results that would be hardly distinguishable from 
the non- retarded curve shown in Fig. 

A further contribution to the vdW force exists from 
the coupling of the field to transverse currents. In fact, 
for good dielectrics and metals and in the region where 
retardation is important, transverse and longitudinal cur- 
rent fluctuations contribute similarly to the vdW force, 
as discussed in section IVl But not for graphene: in ap- 
pendix 1^ this transverse contribution is calculated and 
shown to be, at best, of the order of v/c times smaller 
than the longitudinal part, a result consistent with the 
absence of a retarded regime for the longitudinal part. 

V. SUMMARY: Graphene vs Dielectric and Metal 

Here we summarize graphene's results for the vdW in- 
teraction and, to gain a better perspective, compare them 
with the standard prototypes of metals and dielectrics. 
The basic results for graphene, Eq. [5] and Eq. [TH are 
collected here: 



We have seen that a single length, graphene's thermal 
length: = hv/ksT, controls the vdW force, which ex- 
hibits a crossover from the zero-T results of Dobson et 
al.^ (linked to the linear dispersion of Dirac fermions) 
to a finite-T, universal regime. The latter has been un- 
derstood as arising form the existence of classical plas- 
mons: charge fluctuations of thermally generated carriers 
whose existence and energy scale are tied to tempera- 
ture. These results, originally obtained for the longitudi- 
nal (charge) response with instantaneous Coulomb cou- 
pling, have been shown to survive virtually unaffected 
when retardation is included. In addition, the transverse 
contribution is shown in the appendixElto never compete 
with the longitudinal one. Finally, we have emphasized 
that the temperature dependence of the vdW interac- 
tion in graphene reflects basically a property of matter, 
as opposed to the corresponding thermal regime of the 
archetypal metals and dielectrics that are described in 
the following subsections. 

A. Graphene vs Dielectric 

To make a meaningful comparison, we consider a 
two-dimensional insulator, with characteristic frequency 
(gap) and length scale of atomic dimensions: fi^o ^ eV 
and ao A. We assume low temperatures, so that the 
matter response is well approximated by the zero-T re- 
sponse (room temperature is low temperature for an eV 
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Dielectric Metal 




FIG. 4: Left panel: light-cone {u = cq) versus dynamics of a 
dielectric of typical frequency ujo, illustrating the three vdw 
regimes encountered with increasing distance at finite tem- 
perature (Eq. I21|l . Right panel: As in left panel but for a 
typical 2d metal, with dynamics characterized by plasmons 
(Eq. [22j. 

gap). In this situation, three regimesi^ can be considered 
in the vdW force as a function of separation, correspond- 
ing to the q regions schematically shown in the left panel 
of Fig. [2 

I : / oc —hujQ Oq z < zf^f' 

II: fo,-hcalz-\ 4er^<z<AT (21) 

III : f (X -keTal z-^, z > Xt 

where the characteristic distance for the onset of retarded 
effects is zl'^l'^'^ = c/luq, with the field's thermal length as 
before: At = hc/ksT. Notice that Fig. H and Eq. [H] 
are schematic: transitions between regimes are in the 
form of crossovers, and the very existence of a regime 
(particularly II), and its associated power law, assumes 

well separated values of zf^f' and At- 

Regime I corresponds to the sum of R^^ contri- 
butions expected from a dielectric, where longitudinal 
(charge) fluctuations, instantaneously coupled, dominate 
the force. In region II, matter is quicker than field, and 
retardation matters, implying and additional power with 
distance. Here both longitudinal and transverse currents 
contribute similarly to the vdW force. Finally, a thermal 
regime (III) appears for distances beyond At, at which 
the short-ranged nature of the coupling due the electro- 
magnetic modes with finite Matsubara frequencies leaves 
the mode iujn = as the sole contribution to the vdW 
force. Notice that the T dependence of the vdW force in 
this thermal regime is not due any T-dependence in the 
matter response, which remains that of zero-T for the 
low temperatures considered, as explained before. 

The comparison of Eq. (|20p (and accompanying dis- 
cussion) and Eq. pT|) clearly exhibits that, for graphene, 
it is a difference in the dominant matter response what 
changes the vdW force at finite T for distances above 
the (T-dependent) thermal length, ^t- For shorter dis- 



tances, the zero-T response dominates while at larger 
distances, the coupling of classical charge fluctuations of 
thermally excited carriers dominates. Neither retarda- 
tion nor transverse currents (see appendix E]) do matter 
quantitatively, what explains the absence in graphene of 
a retarded regime like region II of a dielectric. 

B. Graphene vs Metal 

We consider a two-dimensional gas of electrons, treated 
in the RPA approximation and characterized by a Fermi 
wavevector typical of a good metal fc^ ~ A , with Fermi 
velocity of the order of graphene's v ~ c/300 (certainly 
of the order of a typical metal). Interactions are then 
characterized by the same dimensionless number as in 
graphene: a = l(2tohv) ~ 13.6, and just for ease, we 
take iV = 4 species of fermions, as in graphene. As for 
the dielectric, we assume a sufficiently low temperature 
such that the matter response is basically that of zero- 
T (a good degenerate metal). Again, three separation 
regimes can be considere d^"'^^ , with corresponding q re- 
gions schematically depicted in the right panel of Fig. 

a 

/: /oc-Si;fc)/'z-7/2, z<4r*^ 

// : (Casimir) / = - ^ |f , z'^^^^'^ <z<\t 

(22) 

where the characteristic distance for the onset of retarded 
effects in the metal is z^et'^*'' ~ ("^f^^/)"^' ^'^'^ 
field's thermal length as before: At = hc/kBT. Again, 
the regimes (mainly II) of Eq. [221 are meaningful only for 

well separated values of 2:^6*'^*'' ^-^id ^t- 

Regime I corresponds to the instantaneous coupling of 
plasmons, as discussed in refs.[2^ andlsl. It is very impor- 
tant to realize that plasmons, with frequency cjp™'^*'' (q) — 

y^^^v^Itfq, are an intrinsic feature of metals at zero T. 
At the low-T considered here and for wavevectors in re- 
gion I (see Fig. |4]), they are certainly quantum objects, 
;i4""*) » keT. This is in contrast with the situation 
for graphene, where plasmons only exist at finite T and 
with frequency tied to temperature (see Eq. [TT]), there- 
fore they are always classical: hujp^^'^^'^'^"'^^ <C ksT. This 
explains the absence of a vdW regime in graphene, 

whose presence would require free carriers at zero-T, cer- 
tainly not the case in graphene, where this zero-T regime 
is replaced by Dobson et al.^ force. 

Regime II in Eq. [22] is the well known Casimir re- 
sult. Here the field is quicker than matter (plasmons) 
and a fully retarded calculation is needed (see note 
|23| ). with longitudinal and transverse currents in the 
2d metal contributing equally to the Casimir univer- 
sal result. This regime has no equivalent in graphene. 
Again, the reason is the classical behavior of charge 
fiuctuations in graphene. The existence of a Casimir 
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regime would require simultaneously uip 



(graphene) 



(9) » 



cq and fkOff^'^^'^'^"'^^ ^ ksT, something impossible for 
graphene's temperature-tied plasmons (see Eq. [TT|) . The 
quantitative irrelevance of the transverse response for the 
vdW force in graphene, explained in the appendix is 
also consistent with the absence of this Casimir regime 
(transverse currents are responsible of half the Casimir 
result in the 2d metal). 

Finally, thermal effects do appear in regime III, where 
the result for metals is quantitatively identical to that of 
graphene, Eq. [Ml But there are profound differences in 
the reasons that lead to the same expression in both sys- 
tems. In the 2d metal, the matter response remains that 
of zero-T, but temperature forces to select the iuJn = 
term as the only electromagnetic mode that is not short- 
ranged. Therefore, although the result itself, Eq. ([^ 
III), cannot exhibit the dynamics of the field, the explicit 
appearance of the light velocity c in its range of applica- 
bility (z > Xt — hc/ksT) does reveal its origin. In con- 
trast, graphene's finite-T behavior Eq. [T3] (obtained for 
the instantaneous Coulomb coupling) , has been shown to 
arise from the classical plasmonic response of matter, as 
explained before. Therefore the range of applicability of 
the thermal regime in graphene (z > = hv/ksT) does 
not contain the light velocity (infinite for non-retarded 
calculation), but a matter property: Dirac fermion's ve- 
locity (w). 

Let us close mentioning that there has recently been 
much interest in the issue of finite-T vdW interactions in 
poor metals and its relation to non locality in the metal's 
response2ii2^iS&i21. Graphene may well provide a natural 
ground for these concerns as a system exhibiting both a 
zero-T, dispersive-response, result and a classical linear- 
in-T regime, but at much shorter distances than would 
otherwise be required for the classicality of the electro- 
magnetic field: at room temperature, ^ 26 nm for 
graphene versus At ~ 300 for typical metals and di- 
electrics. 
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Vc{q,z) T)t{q,UJ,z) 
fc{q,z) -> ~dz'Dt{q,uj,z) 



(Al) 



where 'Dt{qi^i z) is the (part of the) photon propagator 
that couples to in-plane transverse currents, with expres- 
sion: 



Vtiq.uj, z) 



exp{-q'\z\) 



2 60 q'c 



I „2 



(A2) 



and g' = \J (f' — uP' jf? ■ Xjtjt the transverse current- 
current response of an isolated graphene layer, with RPA 
expression given by: 



Xjtjt w) 



(0) f \ 



l-Vt{q,u;,z = 0)xfX(^,^) 



(A3) 



where the zero-T, non-interacting response is given by 

(A4) 



(0) 2TT 



For graphene parameters, (c ~ 300 u), it is straightfor- 
ward to see that the the RPA treatment is unnecessary 
for the response, therefore we take Xjtjt ~ ^hjt • Trading 
q for g', with the Matsubara sum becoming an integral 
iujn ^ i^^ the resulting expression for the force /t is: 

- / q'^dq' exp(-2<z'|z|) - / dr^- 
27r 7o I" Jo 



(A5) 



where, as in the main text, w^j?^ = v ^ — c ^, and 
a — / {2eohv) . Ignoring the meager renormalization im- 
plied by We//, and to lowest order in v/c, the following 
closed expression for the force is obtained: 



ft 



{NaY vhv 



(A6) 



Although the dependence with distance of Eq. (jA6p is 
the same as that of the longitudinal contribution, Eq. 
m the salient feature is the ratio v/c in the coefficient. 
Indeed, explicit comparison of Ea. (|A6p and Eq. pleads 
to: 



APPENDIX A: TRANSVERSE CURRENTS 
CONTRIBUTION 

Here we present the calculation of the force due to 
the (in-plane) transverse currents, to show that it never 
competes with the longitudinal part. As in section [TVl 
the force per area can be obtained with the following 
correspondence in Eqs. and ([3]): 



/t~ 1.23-/^4.1x10-^/ 

c 



(A7) 



Therefore, the neglect of the transverse contribution 
is fully justified. Even more, our result (Eq. IA7p can 
be seen as a hound for the transverse force over an ex- 
tended range of distances. The reason is that the trans- 
verse contribution here calculated uses the scaling form 
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for graphene's response, Eq. IA41 but this response grows 
unbounded with (imaginary) frequency. Real graphene 
hves in a lattice and, therefore, we expect this response 
to saturate beyond a cut-off of the order of hojo ~ eV 
(the standard diamagnetic local limit). This implies that 
our zero-T calculation ()A6|) should apply for distances 
z > c/wq, while for shorter distances, the saturation of 
matter response must imply even lower (cut-off depen- 
dent) forces from the transverse response. 

Finally, we consider the transverse contribution at fi- 
nite T. It suffices to notice that the zero-T result (|A5[) 
is dominated by the matter response at (imaginary) fre- 
quencies close to the light cone: rj ^ cq. At such high 
frequencies, temperature hardly affects the matter re- 
sponse, and the zero-T transverse contribution holds for 
distances up to z > At (At = hc/kBT). Beyond such 
distances, as usual, only the iun — term survives in 
the frequency sum (Eq. [5]) as the long-range interac- 
tion between graphene's layers. Then, using the long- 
wavelength, finite-T result for the static transverse re- 
sponse: 



xj],(g,iw„ = 0) 



N {v qf 
24^ ksT 



q^T < 1, 



(A8) 



we end up with the following expression for the large- 
distance, transverse contribution to the force per area 
between graphene's layers at finite temperatures: 



ft 



{Na^rfkBT 
f356 7r3 ~^ 



> A7 



(A9) 



This force replaces the zero-T expression (jA6p for dis- 
tances z > At. Quantitatively, it is utterly small when 
compared with the corresponding finite-T longitudinal 
contribution (Eq. I14p . This completes our justification 
for the neglect of the transverse contribution in the body 
of the paper, both at zero and finite temperatures. 



APPENDIX B: REAL GRAPHENE VERSUS 
DIRAC FERMIONS 

The success of the Dirac fermion approach to ex- 
plain the experimentally observed electronic properties 
of graphene leaves little doubt about its correctness for 
describing the low energy physics'^, including room tem- 
perature and above. The long distance behavior of the 
vdW interaction cannot be an exception: qualitative fea- 
tures such as the asymptotic power law of Eq.[8]are linked 
to the linear dispersion of Dirac fermions and its associ- 
ated gapless excitations, as first shown by Dobson et al^. 
On the contrary, sum of treatments would lead to 
an incorrect 1/z^ dielectric law for the force. 



Nevertheless, the presence of further electrons (and 
states) besides those of Dirac cones makes our calcula- 
tion quantitatively incomplete as far as real graphene is 
concerned. These electron-hole excitations being gapped, 
their contribution to the vdW interactions is not ex- 
pected to change the obtained asymptotic behavior. Nev- 
ertheless, estimating the size of this neglected contribu- 
tion is important to assess the range of validity of our 
results. 

Explicit calculation of these contributions is well be- 
yond the scope of this work, but a fair estimation of its 
size can be obtained from published resultsiS,. Such an 
estimation can be restricted to zero T because: i) the 
missing contribution corresponds to gapped excitations 
above ~ 2 eV, hardly affected at any reasonable finite T, 
ii) our results at finite T, Fig. are always greater 

than those at zero T. Therefore, if the missing contribu- 
tions do not compete with the asymptotic regime at zero 
T, they will not at finite T. 

In ref. fisl, an ab initio calculation of the vdW inter- 
action in graphene-like systems is performed. Not in- 
tended for large distances, the low energy physics in 
ref. m is only approximated, leading to an effective long- 
distance interaction between carbon atoms of the tradi- 
tional dielectric form: 1/R~^. Nevertheless, ref. Ts' in- 
cludes electrons beyond the lowest sp band, and high en- 
ergy excitations (up to tens on eV)^'^, providing a quan- 
titative estimate of the missing contributions. The at- 
tractive energy between two atoms in different graphene 
sheets separated by R is obtained there^^ as Ce/R^, with 

A 6 
, and around 50 percent of this con- 
tribution is reported as coming from excitations beyond 
the lowest sp bandit. From these numbers, one obtains 
an estimate of / = —Dq/z^, with Dq ^ 6.3 eV A , 
for the force per area coming from the neglected terms. 
This contribution becomes comparable to that of Eq. [8] 
only at separations between graphene layers of the order 
of 1.5 nm, a rather short distance for a vdW scenario, 
just a few times graphite interlayer distance. Beyond 
this distance, the different power law makes the asymp- 
totic contribution of Eq. [5] dominant (see also ref. |2g). 
For instance, at a distance z ~ 27 nm, corresponding 
to graphene's thermal length a room T, the asymptotic 
contribution here obtained (Fig. [2]) is around 20 times 
the estimate from the ignored terms. At a separation of 
z ^ 100 nm (graphene's thermal length at liquid-nitrogen 
T), the force calculated in this work, would surpass the 
missing contributions ^ 80 times. All this indicates that 
restricting to the low energy physics not only provides 
the large distance asymptotic behavior of the vdW inter- 
action but also the quantitatively dominant contribution 
for real graphene for separation beyond a few nanometers. 
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